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Abstract 

We study the anomalous decays ir°, n — > 77 in the framework of the three- 
flavor Nambu-Jona-Lasinio [NJL] model, in the vacuum and in quark matter in (3 
-equilibrium. It is found that the behavior of the relevant observables essentially 
reflects a manifestation of the partial restoration of chiral symmetry, in non strange 
and strange sectors. The probability of such decays decreases with density, showing 
that anomalous mesonic interactions are significantly affected by the medium. 
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1 Introduction 



The structure of pseudoscalar mesons, its mass spectra and decays have attracted a 
lot of interest along the years, an important motivation for this interest being certainly 
related to the fact that the origin of these mesons is related to the spontaneous or explicit 
breakdown of symmetries of QCD. Furthermore, since at high densities and temperatures 
new phases with restored symmetries are expected to occur, the study of pseudoscalar 
meson observables under those conditions is specially relevant since they could provide 
signs for the phase transitions and the associated restoration of symmetries. As a matter of 
fact, major theoretical and experimental efforts have been dedicated to heavy-ion physics 
looking for signatures of the quark gluon plasma, a state of matter with deconfinement of 
quarks and restoration of symmetries [TJ El EJ • 

In the limit of vanishing quark masses, the QCD Lagrangian has 8 Goldstone bosons, 
associated with the dynamical breaking of chiral symmetry. The non inexistence of a 
ninth Goldstone boson in QCD is explained by assuming that the QCD Lagragian has a 
£7^(1) anomaly. The origin of the physical masses of the different pseudoscalar mesons is 
due to the explicit breaking of symmetries, but it presents differences that will be analyzed 
next. While for pions and kaons it is enough to break explicitly the chiral symmetry by 
giving current masses to the quarks, the breaking of the ?7a(1) symmetry by instantons 
has the effect of giving a mass to rf of about 1 GeV. So the mass of rf has a different 
origin than the masses of the other pseudoscalar mesons and it cannot be regarded as the 
remnant of a Goldstone boson. Concerning the other two neutral mesons, rj and 7r°, they 
are degenerated when the £7a(1) anomaly and the current quark masses are turned off, 
but, when these effects are taken into account it turns out that a percentage of r\ mass is 
due to the anomaly, and, therefore, this meson should regarded as less " Goldstone-like" 
than the pion. Investigating this problem in vacuum and in medium is an important task. 
Besides mass spectrum and meson-quark coupling constants, the observables associated 
with the two photon decay of these mesons might provide useful insight into this problem, 
and calculations of such observables in vacuum and at finite temperature may be found 
in the literature [lEJEJIIllEllIEIiniEEI]. 

Understanding the processes tc°(i]) — ► 77 is specially relevant having in mind that 
the great percentage of photons in the background of heavy-ion collisions is due to the 
decay of 7r°and?7 ^2]- As a matter of fact the production of such mesons is indicated by 
the occurrence of two photon pairs with an invariant mass equal to these meson masses. 
Possible medium modifications of anomalous mesonic interactions is a topic that has 
attracted lost of interest. As pointed by Pisarski [7j and Pisarski and Tytgat jH] while for 
fermions the axial anomaly is not affected by the medium, the opposite is expected for 
anomalous mesonic interactions. In this concern, the study of 7r° — > 77 is particularly 
interesting due to its simplicity and its association with restoration of chiral symmetry. 
Although the lifetime of the neutral pion is much longer than hadronic time scales and it 
is not expected to be observed inside the fireball, the physics is the same as that of other 
anomalous decays (u — > Tmn , uj — ► pir) that are relevant for experiments in the hot/dense 
region. 

A great deal of knowledge on chiral symmetry breaking and restoration, as well as 
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on meson properties, comes from model calculations jTHj. In particular, the Nambu- 
Jona-Lasinio [NJL] [TH ITo] type models have been extensively used over the past years 
to describe low energy features of hadrons and also to investigate restoration of chiral 
symmetry with temperature or density [Tn l lT7|IT ^ IT ^ l271 1 l2T l l2^l2 ^ l2l l l23]. 

This work comes in the sequel of previous studies on the behavior of neutral mesons 
in hot and dense matter (2H1 EH EH] The study of phase transitions in quark matter 
simulating neutron matter in (3 equilibrium, at zero and finite temperature, as well as 
the discussion of behavior of pseudoscalar mesons (in particular neutral mesons) in such 
media, in connection with the restoration of symmetries, has been done in |2*3*1 |2U 122] , by 
analyzing only the mass spectrum. In particular, it was shown that the behavior of the 
masses of 7r° and rj in this concern reflects mainly the restoration of chiral SU(2) symmetry 
and manifests strongly in the pion. This not so evident for the rj, due to its strangeness 
content and the dependence of the anomaly. The aim of this paper is to investigate the 
anomalous decays of 7r° and rj in the vacuum and in quark matter in j3 -equilibrium and 
discuss our results in connection with the breaking and restoration of symmetries. We 
also compare our results with those obtained by other authors who studied the effect of 
temperature on these anomalous mesonic interactions [51 IB1 U\ ITT]. 

2 Formalism 

We consider the three-flavor NJL type model containing scalar-pseudoscalar interac- 
tions and a determinantal term, the 't Hooft interaction, with the following Lagrangian: 

C = q(i^d^ - m)q + \g s £to[UA a g) 2 + {qi^X a qf] 

+ g D {det[g(l + 75) g] + det[g(l - 7b)?]}- W 

Here q = (u, d, s) is the quark field with three flavors, Nf = 3, and three colors, iV c = 3. 
rh = dia.g(m u ,md,m s ) is the current quark masses matrix and A a are the Gell-Mann 
matrices, a = 0, 1, . . . , 8, A = y| I. 

The last term in (JTJ) is the lowest six-quark dimensional operator and it has the 
SUl(3) <£> SUr(3) invariance but breaks the Ua(1) symmetry. This term is a reflection 
of the axial anomaly in QCD and can be put in a form suitable to use the bosonization 
procedure ([201121] and references therein): 

C D = \qd D abc (qX c q)[(qX a q)(qX b q)-3(q il5 X a q)(qz l5 X b q)] (2) 


with constants D abc = d abc if a,b,c G {1,2, ..8}, where d abc are the SU(3) structure 
constants, and D 000 = ^/|, D 0ab = -^J\5 ab . 

The usual procedure to obtain a four quark effective interaction from the six quark 
interaction is to contract one bilinear (qX a q) making a shift (qX a q) — ► (qX a q)+ < qX a q > 
with the vacuum expectation value < qX a q >. Then, an effective Lagrangian can be 
written as: 

£eff = qiij^dft - m)q 

+ S ab [(qX a q)(qX b q)]+ P ab [(qt l5 X a q)(qt l5 X b q)], (3) 
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where 

Sab = gs$ab + gDD abc < q \ C q >, 

Pab = gs^ab ~ gDD abc <q\ c q > . (4) 

Integration over quark fields in the functional integral with (J3J) gives the meson effective 
action 

W eff = -i Tr ln(i(9 M 7 M -m + a a A a + i7 5 a A a ) 

- ^(cr a S ab _1 cr 6 +4> a Pab ~ X 0&)- ( 5 ) 

The fields a a and 4> a are the scalar and pseudoscalar meson nonets. 

The first variation of the action (jSJ) leads to the set of gap equations for constituent 
quark masses Mf 

Mi = m,i -2g s <qiqi> -2g D < q j q j ><q k q k > (6) 

with i ,j ,k = u,d,s cyclic and < >= — iTrSi(p) are the quark condensates. Here the 
symbol Tr means trace in color and Dirac spaces and integration over momentum p with a 
cut-off parameter A to regularize the divergent integrals. The pseudoscalar meson masses 
Mjj (H = Tr , rj) are obtained from the condition (1 - P^W^Pq = M H ,P = 0)) = 0, 
where II 13 (Pq = Mh, P = 0) is the polarization operator at the rest frame 

rF(Po) = 4 [(/j + 4) - (P 2 - (Mi - M,) 2 ) J^(P )] , (7) 
and integrals are given by 

n = £ f (8) 



An 2 Jo E, 

fi(P\ = ^L f A P 2rf P E > + E i | j 1 P* (a) 
2 { 0) 4tt 2 J E % E 3 P 2 - (Ei + Etf + 2ix (E* + E*) ' W 



where Eij = Jp 2 + Mf ■ and E * • = J (p*) 2 + M 2 • are the quark energies. The momentum 



p* is defined by p* = ^(Pq 2 - (Mi - M i ) 2 )(P 2 - (Mi + M i ) 2 )/2P . 
The quark-meson coupling constant is evaluated as 



1 d 

9~Hqq = - Qp (P )]\P =M H , (10) 



where the bound state contains quark flavors i, j. 

Having the on-shell quark-meson coupling constant we can calculate the meson decay 
constant fn according to the definition 

fn = N c g Hqq ^ J ^jtr [(ils)Si(p)( W )Sj(p + P)] . (11) 

Our model parameters, the bare quark masses = m u ,m s , the coupling constants 
and the cutoff in three-momentum space, A, are fitted to the experimental values of 
masses for pseudoscalar mesons (M n o = 135.0 MeV, Mk = 497.7) and f n = 92.4 MeV. 
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Figure 1: The quark triangle diagram for the H — > 77 (direct and exchange process). 

Here we use the following parameterization [2U 1211 125]: A = 602.3 MeV, g s A 2 = 3.67, 
g D A 5 = -12.39, m u = m d = 5.5 MeV and m s = 140.7 MeV. 

We also have M v = 514.8 MeV, 9(M 2 ) = -5.8°, g vm = 2.29, g v - ss = -3.71. 

Note, that 0{M%) is the mixing angle which represents the mixing of A 8 and A com- 
ponents in the 77 -meson state (for details see [2"^ I25j). 

For the quark condensates we have: < uu > = < dd > = — (241.9MeV) 3 , < ss > = 
-(257.7MeV) 3 , and M u = M d = 367.7 MeV, M s = 549.5 MeV, for the constituent quark 
masses. In this section we have described the NJL model with the 't Hooft determinant. 
The model describes well the vacuum properties related with chiral symmetry and its 
spontaneous breaking including their flavor dependence. 

3 The decay H — > 77 

For the description of the decays H — > 77 we consider the triangle diagrams for the 
electromagnetic meson decays. They are shown in Fig. 1. The corresponding invariant 
amplitudes are given by 

f H (P,q 1 ,q 2 ) = ij ■^- i Tr{T H S(p-qi)e 1 S(p)e 2 S(p + q 2 )} 

+exchange. (12) 

Here the trace Tr = tr c trjtr 7 , must be performed over color, flavor and spinor indices. The 
meson vertex function Fh has the ^75 form in the Dirac space, contains the corresponding 
coupling constant guqq (see (I10J1 ) and presents itself the 3x3 matrix form in the flavor 
space. S{p) is the quark propagator S(p) = diag(5' u , Sd, S s ), is the photon polarization 
vector with momentum The trace over flavors leads to different factors for different 
mesons H: Quqq- This factor depends on the electric charges and flavor of quarks into 
the meson H: Q w = 1/3, Q Vu = 5/9 and Q Vs = — y / 2/9. 
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For this evaluation, we move to the meson rest frame and use the kinematics P = qi+q 2 
and P = (Mff, 0). Taking the trace in (|12j) we can obtain 



T H (P, qi , q 2 ) = e^pe^ql T H (P 2 , q\, q 2 2 ) (13) 



where 



and 



%,{P l = M 2 , qf, q z 2 ) = 32a Wflu /^o (14) 



T v (P 2 = M 2 J q 2 1 ,q 2 2 ) = 

3 2o? 7T 

j-^[cos6(5g vilu I^ - 2g vSB fy - smO^hg^uI^ + g^afy]- (15) 



where a is the fine structure constant. The integrals P H = Pjj(P) are given by 

■ f d 4 p 1 

P H (P) = iM t J _ - - _ . (16) 



In order to introduce finite density effects, we apply the Matsubara technique [HJ I2E] , and 
the integrals relevant for our calculation f([5|) and 0) are then modified in a standard 
way [23 • In particular P H {P) (HSJ) takes the form: 

I H (P , P = 0) = -— dp^^-^ ln(— ) (17) 

where Aj is the Fermi momentum. 

Finally, the decay width is obtained from 



p £L T 



and the decay coupling constant is 



gn^rr = % 2!1 - ( 19 ) 



4 Discussion and conclusions 

We present in Table I our results for 7h^ 77 , r#^ 77 and <?h 77 (i? = (7r° ,7/)) in the 
vacuum, in comparison with experimental results (UEIEZI, and we can see that there is 
a good agreement. 

Some comments are in order concerning our calculation of the integral Pjj{P) ()17|) . The 
fermionic action of NJL model (0) has ultraviolet divergences and requires a regularizing 
cutoff. However, the integral P H (P) (fT7j) is not a divergent quantity. We chose to regularize 
the action from the beginning, so that there is no need to cut the integral P H - only the 
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Table 1: Comparison of the experimental values with numerical results obtained in the 
NJL model. 





NJL 


Exp. 




|7;o^ 77 | [eV] 1 


2.5 x 1(T 


ii 


(2.5 ±0.1) x lO" 11 


7T° 




7.65 


7.78(56) 




g^o^n, GcV _1 


0.273 


0.274 ±0.010 






8.71 x 10' 


-IV 


8.57 x 10" iY 






2.54 x 10" 


-li 


(2.5 ±0.06) x 10~ n 


V 




0.440 


0.465 




^-7-7 [GeV] _i 


0.278 


0.260 




T^_» 77 [s] 


1.52 x lO- 


-is 


1.43 x 10~ 18 



ultraviolet integrals I\ (jHJ) and © are regularized [22] • The advantages of using A — > oo 
in non divergent integrals was already shown in |3J 122] , where a better agreement with 
experiment of several observables is obtained with this procedure. 

Now let us discuss our results at finite density. We consider here the case of asymmetric 
quark matter imposing the condition of f3 equilibrium and charge neutrality through the 
following constraints, respectively on the chemical potentials and densities of quarks and 
electrons: p d = p s = p u + p e and \p u - \{p d ± p s ) - p e = 0, with p { = ^(pf - 

Mff' 2 d{n1 - Mf) and p e = 1231231 ES!- 

As discussed by several authors, this version of the NJL model exhibits a first order 
phase transition ^H23I2H]- As shown in [2H], by using a convenient parameterization [21] 
the model may be interpreted as having a mixing phase — droplets of light u , d quarks 
at a critical density p c = 2.25po (where po = 0.17fm -3 is the nuclear matter density) 
surrounded by a non-trivial vacuum - - and, above this density, a quark phase with 
partially restored chiral symmetry [231 12E]- An interesting feature of quark matter in weak 
equilibrium is that at densities above p s ~ 3.8po the mass of the strange quark becomes 
lower than the chemical potential what implies the occurrence of strange quarks in this 
regime and this fact leads to meaningful effects on the behavior of meson observables 

[231 123125]. 

In order to evaluate the transition amplitude of the decay H —>■ 77 in function of 
the density, we require the behavior of Mh and gnqq with density, that are plotted in 
Fig. 2 a)-b), respectively. As we will discuss in the sequel, this behavior is essentially a 
manifestation of the partial restoration of chiral symmetry. The difference between the 
behavior at T = 0,p 7^ and T 7^ , p = is that, in the last case the mesons are 
no more bound states above the critical point since they dissociate in qq pairs at the 
Mott temperature. At finite density they continue to be bound states but with a weaker 
coupling to the quarks as it can be seen in Fig. 2 b). So, it is natural to expect that 
the effect of density on H — > 77 decay observables be qualitatively similar to the effect of 
finite temperature. 

We discuss now our results for the medium effects on the two photon decay of 7r°: 
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Figure 2: Mesonic and quark masses (a) and meson-quark-quark coupling constant (b) 
as function of the baryonic density. 

that are plotted in Fig. 3 a)-c). In order to understand this 
results let us remember the central role of this meson in connection with the breaking 
and restoration of chiral symmetry in the 577(2) sector. A sign for the restoration of this 
symmetry is that the mass of the pion increases with density and this meson becomes 
degenerated with a meson and the pion decay constant f n goes asymptotically to zero. 
The behavior of tt° — > 77 observables with density is closely related with the restoration 
of chiral symmetry in the SU(2) sector. The fact that 7^o_,. 77 , as well as <7„.o_> 77 decrease 
with density reflects the fact that M u and g n Oq q decrease with density. As it can be seen 
from Fig. 2, the quark mass decreases sharply while g n o Uu decreases more slowly. In 
the region p < p c the behavior of the transition amplitude is dictated by a compromise 
between the behavior of the mass and the meson quark coupling. Above p c the mass 
decrease seems to be the dominant effect. Concerning r^o-^ it has a maximum at about 
the critical density, since there are two competitive effects, on one side the decrease of 
the transition amplitude and, on the other side, the increase of the pion mass. Above the 
critical density the two photon decay of the pion becomes less favorable, what reflects that 
it turns a weaker bound qq pair, as already mentioned before. This is compatible with 
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Figure 3: The decay n° — > 77: a) Transition amplitude; b) Decay width; c) Coupling 
constant. 

recent experimental results indicating that pionic degrees of freedom are less relevant at 
high densities |29j . 

In order to clarify the connection between the behavior of the anomalous couplings 
and the restoration of chiral symmetry one can do the simple exercise of calculating the 
transition amplitude in the chiral limit, by setting the external momenta equal to zero. 
Then, similarly to [HE], where such analysis was done with temperature, we get |7^o^ 77 | oc 
^--tt-, the mass decreases being the dominant effect, which leads to a vanishing of the 
transition amplitude at the critical density. 

Concerning the rj — > 77 decay, although qualitatively similar to tt° — > 77, there are 
differences that we will examine in detail and that are related to the evolution of the 
strange quark content of this meson and the behavior of the strange quark in this regime. 
The quark content of 77 is given by: 

1 [2 
1 77 > = cos 9 -j=\uu + dd — 2ss > —sin 9 —\uu + dd+ ss > . (20) 

It was show in j21] that the mixing angle (9 = —5.8° in the vacuum) decreases with density, 
has a minimum (~ —25°) at p ~ 2.8po, equals to zero at p ~ 3.5po ; then it increases rapidly 
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Figure 4: The decay 77 — > 77: a) Transition amplitude; b) Decay width; c) Coupling 
constant. 

up to the value ~ 30°, when strange valence quarks appear in the medium (p ~ 3.8po)- 
g vS s (Fig. 1-b) reflects this evolution of the strange quark content. So, one sees that at 
high densities the 77 is governed by the behavior of the strange quark mass. Since after 
p ~ 3.8p there is a tendency to the restoration of chiral symmetry in the strange sector, 
although less pronounced then for non strange quarks, this effect should manifest itself in 
the behavior of 77 observables. Therefore, although the 7/ is less "Goldstone like" than the 
pion, one expects results qualitatively similar for the two photon decay, which are shown 
in Fig. 4 a)-c). The main difference is that the width T v — > 77 almost vanishes above 
p = 3.8p - This behavior of the 77 seems to indicate that the role of the £/4(l)anomaly for 
the mass of this meson at high densities is less important. 

In conclusion, we studied the behavior of two photon decay observables for the neu- 
tral mesons n° and 77 in quark matter in weak equilibrium and we discussed the results in 
connection with restoration of symmetries. We have shown that, in spite of the different 
quark structure of these mesons, at high densities they share a behavior which is mainly 
a manifestation of restoration of chiral symmetry. We show that these anomalous decays 
are significantly affected by the medium, however the relevance of this results from the ex- 
perimental point of view should be discussed. Recent experimental results from PHENIX 
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j2H] show that 7T° production is suppressed in the central region of Au+Au collisions as 
compared to the peripherical region. This means that tt° — > 77 decay could only be inter- 
esting for experimental heavy-ion collisions at intermediate densities. However, although 
the peak of the tt° width is at a moderate density (see Fig. 3) its life time, r, is here 
of the order of 2.08 x 10~ 17 s, much longer than the expected lifetime of the fireball in 
the hadronic phase, 10~ 22 s, so the decay should occur outside of the fireball. The same 
considerations apply to the 77, although its maximum lifetime is 9.36 x 10~ 19 s. However, 
since the physics underlying these processes is the same of other anomalous processes 
interesting from the experimental point of view, the modification of anomalous mesonic 
interactions by the medium might, in principle, be observed. 
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